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1. Introduction 


Energy is available in multifarious forms and plays a signifi- 
cant role in world wide economic growth and industrialization. 
The growth of world population accompanied with rising mate- 
rial needs intensified the rate of energy usage. Continuous increase 
in energy usage characteristics of the past 50-100 years cannot 
continue indefinitely as demarcated energy resources of earth are 
exploitable. 

On the other hand, environment degradation with use of fossil 
fuels is a menace to life in this earth. In view of world’s deplet- 
ing fossil fuel reserves and environmental threats, development 
of renewable energy sources received importance. Of many alter- 
natives, solar energy stands out a conspicuous energy source for 
meeting the demand. It is considered as adamantine renewable 
energy source due to its huge potential. The freely available solar 
radiation provides an infinite and non polluting reservoir of fuel. 
The easiest way to utilize solar energy for heating applications is to 
convert it into thermal energy by using solar collectors. Solar water 
heaters and solar air heaters are flat plate collectors which used 
mostly for heating of water and air respectively. Solar air heater is 
compact and less complicated in comparison to solar water heater. 
Solar air heater is easy to fabricate with cheaper materials and is 
easy to use then solar water heater. Solar air heater produce hot 
air for any industrial or farmer level drying applications by using 
freely available solar energy, without using any conventional fuels 
like electricity, diesel, LPG, firewood, coal, etc., But it could be cou- 
pled with an existing conventional drying systems like Tray driers, 
tunnel driers, conveyor drier, FBD drier and bin drier operated by 
conventional fuels to save fuel consumption. The thermal efficiency 
of solar air heater is less due to low heat transfer capability between 
absorber plate and fluid flowing in the duct. To make solar air heater 
more efficient solar energy utilization system, thermal efficiency 
needs to be improved by enhancing heat transfer rate. It is reported 
in literature that the heat transfer rate can be exaggerated by the 
following: (a) use of a secondary heat transfer surface, (b) disrup- 
tion of the unenhanced fluid velocity, (c) disruption of the laminar 
sublayer in the turbulent boundary layer. Under the method (c) i.e. 
disruption of laminar sublayer by using wires or ribs in forms of arti- 
ficial roughness provided underside of absorber plate, method (a) 
includes the contact fluid surface area increase by using extended 
surfaces i.e. fins and by using packed bed i.e. a volume of porous 
media obtained by packing particles of selected material into a duct. 
In the packed bed solar energy penetrates to greater depths and 
is absorbed gradually depending on the density of packing. Many 
researchers carried out experimental investigation of heat transfer 
augmentation techniques i.e. extended surfaces, packed bed and 
corrugated absorber in double pass solar air heater but so far no 
work has been reported with artificially roughened double pass 
solar air heater duct. In the present article an attempt has been 
made in order to make a brief outline of the performance double 
pass solar air heater and the various methods those exaggerate their 
performance and mathematical models of some configurations of 
double pass solar air heaters are discussed. 


2. Performance analysis of double pass solar air heater 


It is required to analyze thermal performance of solar air heater 
for making an efficient design of a system. Thermal performance 
concerns with heat transfer process within collector. The different 
arrangement of the double pass solar air is shown in Fig. 1(a)-(c). 
Fig. 1(a) is considered for the thermal analysis and the design 
of such conventional system are described by Satcunanathan and 
Deonarine [1]. 


Nomenclature 

A area element 

B height of solar air collector 
C specific heat 

D depth of air channel 

Dy hydraulic diameter of air channel 
h heat transfer coefficient 

k thermal conductivity 

L length of air collector 

Q Heat collected 

m mass flow rate 

Nu Nusselt number 

Pr Prandlt number 


recycle ratio 
r fraction of mass flow rate 

Re Reynolds number 

S incident solar radiation 

T temperature 

t time 

U loss coefficient 

V wind speed 

W width of absorber 

xy space coordinates 

H fin height 

Z distance from absorber plate 
N number of fins 


Greek letters 

solar absorbtance 
emissivity 

thickness 

reflectance 

transmitivity 

efficiency 

ratio of channel thickness 


D3 adv ere 


Subscript 

ambient 

base plate 
convection 
fluid(air) 
glass cover 
inlet 

outlet 
absorber plate 
radiation 
lower channel 
upper channel 


p =s o mig n D A 


2.1. Thermal performance 


The double pass solar air heater are of mainly two types depend- 
ingon the fluid flow direction namely counter or return flow double 
pass solar air heater and parallel pass double duct solar air heater. 

The analysis of parallel pass arrangement i.e. Fig. 1(c) heat col- 
lected in the lower and upper channel per unit area is defined as 


qı = he, p-(Tp — Ti) + he,1-b(Tp — Ti) (1) 
q2 = hc, g-u(Tp — Tu) + hc,u-p(Tp — Tu) (2) 


where heat collected by lower flow and upper flow are Q; and Qu, 
respectively defined as 


Qi = mCp(Th fo — Tifi) (3) 
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Fig. 1. (a) Schematic diagram of counter flow double pass solar air heater. (b) Schematic diagram of counter flow double pass solar air heater. (c) Schematic diagram of 


parallel pass double duct solar air heater [1]. 


Qu = MuCp(Tu fo — Tufi) (4) 
Total heat transfer Qr = Q; + Q, and efficiency is defined as 

Q 
n= AS (5) 


3. Mathematical models 


A mathematical model of the closed solar air heaters is used 
in particular, to assist in interpreting the observed phenomena in 
the solar air heaters, to design the system, to predict the trends, 
and to assist in optimization. A review of the mathematical models 
for predicting solar air heater systems carried out by Tchinda [47]. 


The models have evolved to a point where several features of the 
process can be predicted, more effort is required before the models 
can be applied to define actual operating conditions as well as to 
further investigate new solar air heater. A review has been done 
on the theoretical modeling of some different designs of solar air 
heaters. Table 1 shows the heat transfer coefficients used in the 
following configured models. 


3.1. Parallel pass solar air heater 


The design of the parallel pass solar air heater is to optimize 
the heat transfer transportation in the solar air heaters. The design 
consists of a glass cover, absorber plate and bottom insulated plate. 
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Table 1 
Heat transfer coefficients used in different mathematical models. 


Reference Solar air heater system 


Heat transfer coefficients 


Forson et al. [15] Parallel pass solar air heater 


Ho et al. [8] Double pass flat plate solar air 


heater with recycle 


Naphon [18] Double pass flat plate solar air 


heater with porous media 


Naphon [28] Double pass flat plate solar air 
heater with longitudinal fins 
Ho et al. [7] Multi-pass flat plate solar air 


heaters with external recycle 
Double pass solar air heater 
with porous media packed 


Ramadan et al. [20] 


El-Sebaii et al. [19] Double pass solar air heater 


with porous media packed 


Radiative heat transfer coefficients from the absorber to the glazing, h, pc and 
from the absorber to the bottom plate, h, pp are given as 

h, pi = (Tp + T?)(Tp + T;)/(€p' +€5' — 1) where i represents c, b. The heat 
transfer by radiation from the glazing to sky is given by 

rcs = 0€-(T2 + T2)(Te + Ts). The forced convection heat transfer coefficient 
between parallel plates is given as: for laminar flow (Re < 2300) [36] 

(HD, /k) =5.4 + 0.00190(RePrD,,/L)!-7!/1 + 0.00563(RePrD,,/L)'-7!, for transition 
flow region (2300 < Re < 6000) [37]: 

(HD, /k) =0.116(Re2/3 — 125)Pr'/3[1 + (Dp/L)27 (u/uw)®4 and for turbulent flow 
(Re > 6000) [38] (hD;,/k)=0.018Re° Pr, hrs = o&g(Te — T$ )/(Tg — Ts). Free 
convective heat transfer coefficient is given by: 

h = k/Dp,(0.0158Re?,* + (0.0181Rep, + 2.92) exp(—0.0379x/D},) where 

Rep, = 2M/(1 + D/W). Up =[1 + 2(4in +D + 5g/W)][5in/Kin + 1/hw]-1. 

Ug = (UA) eage [Ac = {(k/x)S;}/W. 


hy = h}, hy = hy. The forced convection heat transfer between two plates see 
Kays [39], for short conduit hı and hz see [40]: 

Nu; = hjDe /k = 0.0158Re?*[1 + (De /L)™ |G = 1, 2), for the laminar flow see 
Heaton et al. [36], 

1/hy = 1/heis + 1/h, + lg/kg 1/Ugis = 1/hrgzs + hw + 1/Hegig2 + hrgig2 where 
Negiga a 1.25(Tgım = Tg2m pa see [41 ] 

hpg = 4oT? /[A/ep) + (1/er) —1], hrr = 4073, /[/ep) + (1/er) — 1], those 
between two glass covers see [42] and from cover 2 to ambient, 


hrg2s = o£g(T2, + Tiy XTg2 + Tsky). 


The radiative heat transfer coefficients see [15]. The convective heat transfer 
coefficient for air flowing over the outside surface of the glass cover is 
proposed by [40] hy =5.7 +3.8 V, convective heat transfer coefficient for 
flowing in the channel is given by [43]: Nu=(hD,/k)=0.333Re°8Pr!/3 with 
Dp =4WB/2W + 2B. 


The convective heat transfer coefficient between channel see [28], for radiative 
heat transfer coefficient see [15]. 


hy =h}, h2 =h}, hg =h} and hy = h}. Other coefficients see Ho et al. [8]. 


Hegipt NUpiugi /De = 0.2Re0 Pr"? ky, /De see [44], with Rep =myDe/ Api, 

De =2/3[€Dpi/(1 — €)], £ is the porosity of packed bed (e <1): 

€=(Vohan — Vp1 )/Vchan, Vehan and Vp; are the total volume of the channel and the 
volume of the packed bed material. Dpı = 6VP nx see [45], Vs = total volume 
of n particles selected randomly. Ap =6(1 — €)/Dpi see [45], hpr =Ncgif: see [15], 
hepıfi =[1/hpıf + De/Dp1Sp1 J see [46] with Npıf = hpıfDelkfi = (0.255/2) 

Re? pr", Spi is the constant depends on the packed bed shape. 


The heat transfer coefficients see [20]. 


There are two air flow channels which one is located between the 
glass cover-absorber and another is located between the absorber 
and the bottom insulated plate as shown in Fig. 1(c). Forson et al. 
[15], obtained the energy transport equations for the glass cover, 
channel 1, absorber plate, channel 2, and the base plate is following: 


At the cover: 


oT. 
MgCg FE = gS + hrpg(Tp — Tg) + hef-1g(Tf1 — Tg) 


hcgw(Tg Tw) hrga(Tg Ta) (6) 
For air flowing between cover and plate absorber: 


dT GC OT 
Mpap at W ax 
For plate 1: 


t Nepp1 (Tp — Tei) — Mefig (Th — Tg) (7) 


aT, PT, 
Mi Cpr = OpteS — kpôp > > — Ripe Tp — Te) — hepp2(Tp — Tha) 


—hypp(Tp —Tp)- Nepp1 (Tp — Tj) (8) 
For air flowing between plate absorber and base plate: 


Ty.  —G2Cp ðTp2 
Mpap ot W ox 


+ hepf2(Tp — Tf2) — heof2(Tp — Tr2) (9) 


For bottom plate: 


oT, 07T, 
MpCp a = kp dp aE hrpb(Tp Tp) 


hepf2(To — Tf2) — hp(Tp — Tr) (10) 


Along with following boundary and initial conditions 


Tp 
Ox 


a 
z9: Ox 


x=0 x=L 


3.2. Counter flow solar air heater 


Ho et al. [8] developed the mathematical model for the counter 
flow solar air heater is shown in Fig. 2. In such collectors air is 
allowed to move faster than the simple collectors. The following 
assumptions were considered. 


a) Temperatures of absorbing plate, bottom insulated plate and the 
fluid are only functions of air flow direction. 

b) Glass cover and absorber plate do not absorb radiant energy. 

c) The radiant energy absorb by the outlet cover is negligible. 
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Fig. 2. Multipass flat plate solar air heater with recycle [7]. 
Considering above assumptions, energy equations of cover 1, [4(1 — A)H]w 


absorber plate, bottom plate, channels are as following: 


For cover 1: 
hrpg2(Tp — Tg1) + A (Tp(Z) — Tg1) = Ueis(Te1 — Ts) (11) 
Absorbing plate: 
OpTgiTg25 — hr(Tp — Ts) + hg(Tp — Ts) + hı (Tp — Ty (x)) 

+ (Tp — Talx) i 
Bottom plate: 
Arpr(Tp — Tr) + h4 (Ta(x) — Tr) = Upis(Tp — Ts) 13) 
Bottom channel: 
ha(Tp — Ta(X)) — hy(Ta(x) — Ta) = Ë E Te 14) 
Upper channel: 
hy (Tp — Tp(x)) — h (T(x) — Tg1) = (R+1)mCp dTp(x) (15) 


W dx 


In this the characteristic length is the equivalent diameter of the 
duct. For open conduit, lower channel, upper channel, they are 
respectively given by: 

4AHW 
2AH+W 


Dea = 


and Dep = (16) 


[2(1 — A)H]+ W 
and the air velocities of the open conduit, lower channel, upper 
channel used by Ho et al. [8] is given by: 

(1+R)m j [m(1 +R)] 
= 2AHWp’ ”~ RO- AH]Wp 


m 
= 2HWp’ 


Vo Va (17) 


3.3. Double pass flat plate solar air heater with porous media 


Naphon [18] studied theoretically the heat transfer characteris- 
tics and performance of double pass flat solar air heater with and 
without porous media. The mathematical model proposed in this 
study of design of the double pass single solar air collector with 
porous media is based on the model proposed by Naphon and Kon- 
grtagool [6] and Wijeysundera [2]. The energy equations of the 
system are as follows: 


At cover glazing: 


AgS = hrgs(Tg — Ts) + heg-f1(Tg — Tf1) + hw(Tg — Ta) 


+ hrgp(Tg — Tp) (18) 
At absorber plate: 
QpTgS + hrgp(Tg — Tp) + hefip(Tfi — Tp) = hepf2(Tp — Tp2) 

+ Hrpp(Tp — Tp) (19) 


486 S. Chamoli et al. / Renewable and Sustainable Energy Reviews 16 (2012) 481-492 


At rear plate: 


hef2b(Tf2 — Tp) + Arpp(Tp — Th) = Up(Tp — Ta) (20) 
At a fluid between absorber plate and a rear plate: 

aTpo d? Tro 
mO -i = K- + hef2p(Tp — Tf2) + hef2b(To — Tf2) (21) 


The implicit method of finite difference scheme had been employed 
to solve these models. 


3.4, Double pass flat plate solar air heater with longitudinal fins 


Naphon [28] has investigated heat transfer characteristics and 
entropy generation of the double pass flat plate solar air heater 
with longitudinal fins as shown in Fig. 4. The mathematical model 
developed with the energy balance equations as follows: 


For glass cover: 


gS = hw(Tg — Ta) + hefice(Tg 


T1) + hrgp(Tg — Tp) + hrag(Tg — Ta) 


(22) 
For the first air pass: 
mG -ix = hgig(Tp — T1) + hef1p(Tp — Tri) 
N Z=H 
l 2Lhı(Tyı — Tp) dz (23) 
Afrontal Z-0 
For upper fin: 


Z=H 
N ( am ) N | 
kA = 2Lhy(Ty1 — Try) dz 24 
Afrontal dz /z~0 Afrontal 7-0 (Tn n) (24) 


For absorber plate: 


Og TgS = hef1plTp — Tf1) + hef2p(Tp — Tr2) + hrgp(Tp — Tg) 


N dT, dT, 
Tp) 4 ( kA 1) j ( kA 2) 25 
b) Afrontal dz /z-0 sf dz es) 


t Aypp(Tp 


For lower fin: 


Z=H 
N ( ay N | 
kA = 2Lhz(Ty2 —Ty2)dz (26 
Afrontal dz }z-0 frontal Jz- 2(Ty2 — Ty2) (26) 


For second air pass: 


dTi 
mG- = hgaglTp — T2) + hefaplTp — Tha) 
N Z=H 
L 2Lh2(T,2 — Ty2) dz n 
Afrontal J0 


For bottom plate: 


Tr2) t h,pp(Tp 


hef2p(Tp Tp) + Ua(Tp — Ta) = 0 (28) 


3.5. Multi pass flat plate solar air heaters with external recycle 


Ho et al. [7] have developed a theoretical formulation for a multi 
pass solar air heater as shown in Fig. 2 with external recycle and 
investigate the recycle effect on efficiency. In this study Ho el al. 
have considered a recycled four pass solar air heater with verti- 
cal thickness 2H and 2(1 — œ)H (œ is the channel thickness ratio 
in vertical direction) and horizontal width 26H and 2(1 — £)H (6 is 
the channel thickness ratio in horizontal direction). The following 
assumptions were considered for developing mathematical model: 
(a) the temperatures of the absorbing plate, bottom plate and bulk 


fluid were functions of the flow direction only, (b) both glass cov- 
ers and fluid does not absorb radiant energy, (c) the glass cover 2 
absorbtance is negligible, the energy equations are as follows: 


Cover1: 
hrpgı(Tp — Tg) + h3(Te(x) — Tg1) + h4(Ta(x) — Tg1) = Ugis(Tg1 — Ts) 
(29) 
Absorbing plate: 
QpTe1 Tg2S — hr(Tp — Ts) — hg(Tp — Ts) — hı (Tp — Ta(x)) 
— h2(Tp — Ty(x)) + h3(Te(x) — Tp) + h4(Ta(x) — Tp) = 0 (30) 


Bottom plate: 
hı(Ta(x) — Tr) + h3(Tp(x) 


Tr) +4 hrpr(Tp Tr) = Ugs(Tp Ts) (31) 


Channel a (flowing under the absorber plate) : hı(Tp — Ta(x)) 


ET Ere] ae (32) 
Channel b : h2(Tp — Tp(x)) — h' (Tp(x) — Tr) 
se] : 
Channel c : h3(Tp — Te(x)) — h4 (Te(x) — Tg1) = [o] ka 
(34) 


Channel d : h4(Tp — Ta(x)) — h4(Ta(x) — Tg1) = | RmCp | dTa(x) 


Ww(1-£ß)] dx 
(35) 
In this study, equivalent diameters are given as: 
4aHWB 4aHW(1 — £) 

ea = aH +WB’ ©? 7 2aH+W( — By’ 
D 4(1 — a)HWB 4(1 — a)HW(1 — £) (36) 

ec = 20 —-a@H+ wp’ 17 21-a) +W- B) 
and average velocities are given as: 
v ma ri Mp v ne 

= 2aHWBp’ >” 2aHW(1— pjo’ ~ 21 —a)HWBp’ 
va = ma (37) 


2(1 — a@)HW(1 — B)p 


The above equations were solved numerically using an iterative 
procedure. 


3.6. Double pass solar air heater with packed bed 


The schematic diagram of double pass solar air heater with 
packed bed in the upper channel is shown in Fig. 3 Ramadan et al. 
[20] investigated heat transfer in double pass solar air heater with 
porous media in upper channel. The mathematical model has been 
developed taking into account the following assumptions: (a) the 
air heater operates under steady state condition, (b) the heat capac- 
ities of the glass covers, absorber, and back plates and insulation are 
negligible, (c) the temperature of flowing air was varied only in the 
flow direction. By taking account of all the above assumptions the 
energy balance equation written as follows: 
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“hy == -Y N honin Ssasa-r= Cover 1 
ne hp o LE nan hepri Absorbing plate 
hepr2 hpb =) 
hepr2 Bottom plate 


Fig. 3. Schematic view of energy flow of double pass with packed bed in upper channel [20]. 


Cover 2: 
AgægS = Ag(hrgig2 + Negigi2)(Tg2 — Tg1) + hwAg(Tg2 — Ta) 
+Aghygosky(Tg2 — Tsky) (38) 
Cover 1: 
Ag@gTgS = Ag(hrgig2 + hegig12)(Tg1 — Tg2) — MypigiAg(Tp1 — Tz1) 
— hegıf1Ag(Tf1 — Tg1) (39) 


Air flow in upper channel: 


dTa 
Whpf1(Tp — Tfi) + Whepifi(Tpi — Tri) = mG- 
+ Whegıf1(Tfı — Tg1) + WUs(Ty1 — Ta) (40) 


Absorber plate: 


Aphrp1p(Tp1 — Tp) = Aphpfı(Tp — Tja) 
+Aphpf2(Tp — Tf2) + Aph,po(Tp — Tp) (41) 


S 


f hepifi 


— 


Air flow in lower channel: 


ars 
Tyo) = Mf2Cf2 Fi + WUs(Tf2 — Ta) 
(42) 


Whpf2(Tp — Tr2) + Whf2b(Tp 


Back plate: 
Abhpf2(Tf2 — Ty) + Abhrpp(Tp — Tp) = ApUp(Tp — Ta) (43) 


The various heat transfer coefficient given in the above equations 
were calculated using correlation given by [42]. 

El-Sebaii et al. [19] have proposed an investigation of heat trans- 
fer on a double pass solar air heater in which air is first flow through 
upper channel and then circulated in lower channel packed with 
porous media as shown in Fig. 4.The energy balance for the above 
system were as follows: 


Cover 2: 
AgdgS = Ag(hrg1g2 + hcg1g12)(Te2 = Tg1) F hwAg(Tg2 — Ta) 
+Aghygosky(Tg2 — Tsky) (44) 
t hyg2sky 
Cover 2 
hyot 2 
Cover 1 
hegi f1 


Absorber Plate 


Ppip hepe2 J l hrpb 


EEE 


Bottom Plate 


hepe2 


Fig. 4. Schematic view of counter flow solar air heater with porous media in lower channel [19]. 
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Cover 1: 
AgtgTgS = Ag(lrgiga + Negig2 XTg1 = Tg2) = —hrpıg1Ag(Tp1 = Tg1) 
—Aghoi fi (Tfi — Tg1) (45) 


Air flow in upper channel: 


Whop1 (Tp — Tfi) + Whepifi(Tpi — Tri) = Mpa ot 
+ Whegipi(Tp1 — Tg1) + WUs(Tp1 — Ta) (46) 
Absorber plate: 
Op Tg”ApS = Aphepp1 (Tp — Tr) + Aphepra(Tp — Tra) 
+ Aplypet (Tp — Tg1) + Aphrpp1 (Tp — Tpi) (47) 


Air flow in lower channel: 


Whpf2(Tp — Tro) + Whepif2(Tp1 — Tf2) = Mf2Cf2 oe 

+ WUs(Tp2 — Ta) + Whpf2(Tf2 — Tp) (48) 
Back plate: 
Aphpf2(Tf2 — Th) + Aphrpip(Tp1 — Tp) = ApUp(Tp — Ta) (49) 


The various heat transfer coefficients in the above equations 
were calculated using the correlations given by [44,45]. 


4. Smooth double pass solar air heater 
4.1. Counter flow solar air heater 


Satcunanathan and Deonarine [1] introduced the idea of dou- 
ble pass counter flow air heater. The air flows first between the 
space between glass cover and absorber plate and then induced 
through the duct. The thermal losses from one or more cover glass 
can suppress by using such systems. These systems are 10-15% 
more efficient then single pass system. A heat transfer model for 
two pass system is given by Wijeysundera et al. [2] and results 
were compared with single pass system. It is encountered that the 
open systems, with inlet fluid at ambient temperature, the dou- 
ble exposure systems are 10-15% more efficient then the single 
exposure for a wide range of operating conditions. It is articulated 
that the outer cover temperature is conterminous to ambient tem- 
perature compared to the system operated on single cover glass. 
A mathematical model for performance prediction of a two glass 
cover operated under single and double pass mode developed by 
Persad and Satcunanathan [3]. Study reveals that the performance 
of double pass is the independent of the length, over collector range 
considered and the critical plate spacing depends on temperature 
of the glass cover of collector. Chaudhary and Garg [4] modified 
the mathematical model proposed by Wijeysundera et al. [2] in 


Glass cover Absorber plate 


hÅ... nn SS SSS ——————— es 
= 2(1-A)H_ (1+R)m b channel 
SS ESN) 
Ta 2AH (1+R)m a channel 


order to study the effects of parameters like duct length, upper and 
lower air flow duct depths and mass flow rate of fluid in the sys- 
tem. It is found that the optimum duct depth, corresponding to 
minimal annual cost per unit solar energy gain, were different for 
different duct lengths and air flow rates. Recycling concept in air 
heater is introduced by Ho and Yeh [5] and reported that increasing 
the fluid velocity by using recycle in double pass parallel plate air 
heater enhance the heat transfer coefficient, resulting in improved 
performance. On increasing the fluid velocity, recycling produces 
an effect of remixing the inlet fluid with hot outgoing fluid. The 
mathematical model of Wijeysundera et al. [2] used by Naphon 
and Kongtragool [6] to field out numerically the effect of air flow 
rate on heat transfer characteristic and performance of the single 
and double pass solar air heater. The explicit technique of FDM was 
used to solve the models. Performance of multi pass solar air heater 
with external under countercurrent flow operations was studied by 
Ho et al. [7]. Study shows that the recycle effect increases the col- 
lector coefficient due to the increasing fluid velocity effect and may 
compensate for the temperature difference decrement. The recycle 
effect increases convective heat transfer rate and reduce the tem- 
perature driving force and thus the thermal performance increases 
and further improvement in collector efficiency is obtained when 
channel thickness ratio deviates from 0.5. The recycle effect on dou- 
ble pass flat solar air heater also studied by Ho et al. [8] is shown in 
Fig. 5 found an increase in convective heat transfer rate and reduces 
the temperature driving force and thus increases thermal perfor- 
mance. It explained the two effects originated with recycle; the first 
associated with the fluid mixing of hot and cold air which increases 
heat loss and the second is high rate of heat transfer thus the effect 
overcome the come and tends to enhance thermal performance 
of solar air heater, the performance decreases when the channel 
height deviates from 0.5 and 28-95% increase in collector efficiency 
improvement in comparison to the sub collectors in series. 

A theoretical investigation of double pass solar air heater with 
constant heat flux by using an eigen function expansion in terms 
of power series for the homogeneous part and an asymptotic solu- 
tion for the non homogeneous part was carried out by Ho et al. 
[9]. Study shows that the heat transfer efficiency increases with 
increasing graetz number and as sheet location is moved from 
0.5. The mathematical model can be used in any geometry with 
more general boundary conditions of multi stream or multiphase 
problem and to predict the wall temperature distribution of the 
double pass systems. Ho et al. [10] experimentally and theoreti- 
cally investigates a system where absorber plate divides the duct 
into two channels provided with external recycling. Experimental 
and theoretical results are represented graphically and compared 
with that in the downward-type single-pass solar air heaters of 
the same size without recycling. Considerable improvement in heat 
transfer was observed by employing double-pass operations with 
external recycling and fin attached over and under the absorbing 
plate. The influences of recycle ratio and absorbing plate location on 


Base plate 


Fig. 5. The recycle effect on double pass flat solar air heater [8]. 
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the heat-transfer efficiency enhancement as well as on the power 
consumption increment have been also delineated 


4.2. Parallel pass double duct solar air heater 


In parallel pass solar air heater the air flow is from upper and 
lower of the absorber plate simultaneously, in which there are 
reduced thermal losses in comparison to conventional system. The 
non porous double flow solar air heater system studied by Sodha 
et al. [11] and the study reveals that when air flows above and under 
the absorber plate, thermal performance is better comparative of 
single pass system. Mathematical model and solution procedure 
for predicting the thermal performance of a solar air heater for four 
different configurations: Type-1 single channel with single air blow 
between top glass and bottom absorber plate, Type-2 single chan- 
nel with single air flow between absorber and base plate with no 
glass cover, Type-3 double channel with single air flow between 
absorber and bottom plate and Type-4 double channel with double 
air flows between glass and absorber plates was developed by Ong 
[12]. The model predicts the mean temperature of the heated air 
near the inlet of the collector and under predicts near the outlet. An 
experimental investigation in order to predict the thermal perfor- 
mance of a system with simultaneous air flows over and under the 
absorber plate was carried out by Yeh et al. [13], in order to study 
the fraction of mass flow rate in upper and lower channel. A consid- 
erable improvement in collector performance of double pass over 
single pass at same flow rate and in order to maximize the thermal 
performance of double pass collector the mass flow rates in both 
channel must be same. A model of single pass double duct solar air 
heater is given by Hegazy [14]. The study reveals that the channel 
depth to length ratio should be 0.0025 for both free and forced con- 
vection. The optimum value is considered and a theoretical study 
on collector shapes is conducted and the results shows for a con- 
stant area, variables width collectors exhibit performance behavior 
similar to constant width ones, but with a marginal decrease in col- 
lector efficiency. Forson et al. [15] Developed mathematical model 
for single pass double duct solar air heater (SPDDSAH) which will 
be a design tool capable of predicting the incident solar radiation, 
heat transfer coefficient, mean air flow rates, mean air temperature 
and relative humidity at the exit. It was found that on increasing 
the top to bottom channel ratio is a cost effective way of optimizing 
the thermal performance of SPDDSAH. 


5. Double pass packed bed systems 
5.1. Double pass solar air heater with packing material 
There are two problems associated with the thermal perfor- 


mance of solar air heater. One is the heat transfer from the metallic 
plate to the flowing air is very less as low thermal conductivity of 


Absorber 


air. The second associated with the low heat transfer due to small 
thickness of absorber plate. This is increased by using a porous 
bed absorber. As the voids exits in the porous bed, solar radia- 
tion impinges to a greater depth and thus relatively absorbs there. 
The packing provides an increase in the turbulence which results a 
high rate of heat transfer. The packing material also acts as a heat 
storing media. A numerical study given by Mohamad [16] for inves- 
tigation of the performance of double pass counter flow solar air 
heater with and without porous packed bed and its comparison 
with single and double glazing conventional solar air heaters. It 
was reported that thermal efficiency of a counter flow air heater 
without porous matrix is between 10 and 18% higher than the dou- 
ble glazing conventional air heater. The thermal efficiency of single 
glazing conventional air heater is between 18 and 25% lower than 
the counter flow air heater without porous matrix. An indoor test 
conducted by Sopian et al. [17] on the double pass solar collec- 
tor with and without steel wool as porous media in the second 
channel. The collector had only single glass cover and a black- 
ened metal absorber and the material used as the porous media 
was steel wool. Various combinations of upper and lower channel 
were considered for the investigation. The porous media in the sec- 
ond channel results an increase in outlet temperature and 20-60% 
increase in the thermal efficiency was reported. The heat trans- 
fer characteristics and performance of double pass solar air heater 
with and without porous media in the lower channel was studied 
by Naphon [18] as shown in the Fig. 6. The mathematical model 
developed in the study based on the systems of Wijeysundera et al. 
[2]. The effect of thermal conductivity of the porous media on the 
heat transfer characteristics was studied and the results validated 
with experimental results of Sopian et al. [17]. 

An experimental investigation of the thermal performance of 
double glass, double pass solar air heater with packed bed (DGDP- 
SAHPB) was carried by ElI-Sebaii et al. [19] to study the effect of 
mass flow rate and porosity on the air outlet temperature, ther- 
mal output power, pressure drop and thermo hydraulic efficiency. 
Study reveals that the heat transfer is enhanced by using gravel 
instead of lime stone as a packed bed above the absorber plate 
and the annual averages of the outlet temperature and thermo 
hydraulic efficiency of the DGDPSAHPB are about 16.5% and 28.5% 
higher than those of DGDPSAH. A theoretical and experimental 
investigation of double glass double pass solar heater with packed 
bed (DPSAHPB) above the heater absorber plate carried by Ramadan 
et al. [20]. The lime stone and gravels were used as packed bed 
materials and effect of mass flow rate of air and porosity of packed 
bed material was studied and it was found that for increasing the 
outlet temperature of the flowing air after sunset, it is advisable 
to use higher masses and low porosities. The thermo-hydraulic 
efficiency was found to increase with increasing mass flow rate 
until the flow rate of 0.05 kg/s beyond which the increase in effi- 
ciency becomes insignificant and operating range of mass flow rate 
is set to 0.05 kg/s or less in order to have low pressure drop across 
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Fig. 6. Double pass solar air heater with porous media in the lower [18]. 
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Fig. 7. Double pass solar air heater with longitudinal fins [28]. 


the system. Theoretical and experimental investigation of double 
pass solar heater with porous material packed in the lower channel 
carried by Sopian et al. [21]. The porous media arranged in differ- 
ent porosities to increase heat transfer, area density and the total 
heat transfer rate. Heat transfer and pressure drop relations have 
been developed for air flow through porous media. Porous media 
in the second channel increases the heat transfer coefficient and 
thus improved thermal efficiency. The thermal performances of sin- 
gle and double pass solar air heaters with steel wire mesh layers 
instead of a flat absorber plate were investigated experimentally 
by Aldabbagh et al. [22]. The results indicate that the efficiency 
increases with increasing the mass flow rate with in the range 
of 0.012-0.038 kg/s. For the same range of flow rate, efficiency of 
double pass is 34-45% higher than single pass and the maximum 
efficiency attained in single and double pass is 45.93% and 83.65%, 
respectively for the flow rate of 0.038 kg/s. Thermal performance of 
single and double pass solar air heater with fins attached and steel 
wire mesh as absorber was carried out experimentally by Oma- 
jaro and Aldabbagh [23]. Effect of air mass flow rate over the range 
of 0.012-0.038 kg/s on the outlet temperature and thermal effi- 
ciency was studied for the bed height of 7cm and 3cm of lower 
and upper channel, respectively. The efficiency of double pass was 
7-19.4% more than of single pass system. The maximum efficiency 
attained at the flow rate of 0.038 kg/s and found 59.62% and 63.74% 
for single and double pass system, respectively. The proposed sys- 
tem of Mohamad [16], Aldabbagh et al. [22] is modified by Ramani 
et al. [24] and studies carried out on the system with air forced to 
flow in the lower channel packed with wire mesh matrix as porous 
media. The predicted values of the system validated with the out- 
door studies and it was found that the efficiency of double pass 
with porous is 25% higher than without porous and 35% higher than 
the single pass system. An experimental study of double pass solar 
air heater with packing material on the upper channel was car- 
ried out by Dhiman et al. [25] and developed an analytical model 
was developed that describes the various temperatures and heat 
transfer characteristics of parallel flow packed bed solar air heater, 
employed to study the effects of mass flow rate and varying porosi- 
ties of packed bed material on its thermal performance. PFPBSAH 
thermal efficiency was 10-20% higher than the conventional dou- 
ble pass system. The thermal performance of a double pass solar air 
heater with 2, 4, and 6 fins attached was investigated experimen- 
tally by El-khawajah et al. [26] and wire mesh layers were used 
between the fins instead of an absorber plate. The results indi- 
cate that the efficiency increases with increase in mass flow rate 
for the range of the flow rate used in this work between 0.0121 
and 0.042 kg/s. Moreover, the maximum efficiency was obtained 
by using 6 fins at the same mass flow rate. The maximum effi- 
ciency obtained for the 2, 4, 6 fins of SAH were 75.0%, 82.1% and 
85.9%, respectively for the mass flow rate of 0.042 kg/s. In addition, 
the maximum average temperature difference between the inlet 
and the outlet, AT, for the SAH with 6 fins was the highest for the 
same mass flow rates compared to 2 and 4 fins SAH. The maximum 
average and instantaneous peak of AT obtained were 43.1°C and 


62.1 °C, respectively for the 6 fins SAH when the mass flow rate was 
0.0121 kg/s. 


6. Double pass system with extended surfaces 
6.1. Counter flow with extended surfaces 


Experimental investigation of heat transfer characteristics of a 
double flow channel with fins attached on the upper and lower sides 
of absorber plate was carried out by Yeh et al. [27] and fins of rectan- 
gular profile of height 5.5 cm and pitch of 6 cm were attached to the 
absorber plate. Both the theoretical and experimental investigation 
show that the optimal ratio of air flow rates in both the flow chan- 
nels for maximum efficiency is unity i.e. both the channels have 
same mass flow rates. The double pass air heater with fins attached 
show a considerable increase in the performance as compared to 
the smooth one. A numerical study on the performance and entropy 
generation of the double pass solar air heater with longitudinal fins 
shown in Fig. 7 was carried by Naphon [28] and a mathematical 
model was developed for the prediction of heat transfer charac- 
teristics for the range of mass flow rate 0.02-0.1 kg/s. The effect of 
height and number of fins on entropy generation was considered 
and it was found that thermal efficiency increases with increasing 
the height and number of fins whereas the entropy generation was 
inversely proportional to the height and the number of fins. 

Ho etal. [29] Carried out experimental and theoretical investiga- 
tion of a double pass solar air heater with fins attached as shown in 
Fig. 8 and the experimental and theoretical results were compared 
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Fig. 8. Double pass solar air heater with fins and baffles [29]. 
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Fig. 9. Solar air heater with internal fins attached [30]. 


with downward type single pass solar air heaters of the same size 
without recycling. The study reveals that the baffled double pass 
solar air heater is a feasible design to improve collector efficien- 
cies with recycles operations. As baffle create high turbulence and 
extending the heat transfer area and thus a high rate of heat transfer 
coefficient is achieved. It was noticed that on increasing the mass 
flow rate and recycle ratio is found to be beneficial for collector effi- 
ciency improvement; higher energy dissipation occurs at the same 
time when operating at high mass flow rate and recycling ratio. 

A theoretical investigation has been carried out by Yeh and Ho 
[30] in order to study the influence of external recycle on the effi- 
ciency of solar air heater with internal fins attached to it is shown in 
Fig. 9. The external recycle produces two effects one is the increase 
in fluid velocity to decrease the heat transfer resistance while the 
second one is the lowering the driving force of heat transfer, due 
to remixing at the inlet, that a bad for performance. It was found 
that there will be considerable amount of efficiency increment if 
the operation is carried out with external recycle to suppress the 
undesirable effect. It was noted that enhancement increases with 
increase in reflux ratio, especially for operating at low flow rates 
with higher inlet temperatures. A theoretical investigation of the 
collector efficiency of upward type double pass solar air heater with 
fins attached and with external recycling was carried by Ho et al. 
[31] and a mathematical model was developed for predicting the 
heat transfer characteristics of the system. The equations for the- 
oretical prediction of the useful heat gain, collector efficiency and 
outlet temperature were derived. It is found that desirable effect of 
increasing the fluid velocity using recycling to overcome the unde- 
sirable driving force decreases (temperature difference) for heat 
transfer due to remixing at the inlet, more than 80% improvement 
in collector efficiency is obtained using recycling operation. 


6.2. Parallel pass flow with extended surfaces 


Experimental energy and exergy analysis for a novel flat plate 
solar air heater with different obstacle in fluid flow was carried 
out by Esen [32], the operating parameters that are to be mea- 
sured for double flow air heater were inlet and outlet temperature, 
the absorbing plate temperature, the ambient temperature and 
the solar radiation. It was found that double flow solar air heater 
increases the heat transfer area and thus the heat transfer rate, 
thermal and exergetic efficiencies. 

Ozgen et al. [33] experimentally investigated absorber plate of 
double pass channel solar air heater. Three arrangements were used 
for absorber plate: (a) cans with staggered arrangement, (b) cans 
with inline arrangement and (c) plane duct without cans is shown 
in Fig. 10. The range of mass flow rate over which experiments per- 
formed was 0.03 and 0.05 kg/s. The double flow solar air heater with 
aluminum cans increases the heat transfer area and thus the con- 
vective heat transfer and the optimum value of collector efficiency 
was obtained for staggered arrangement at the mass flow rate of 
0.05 kg/s. 

Theoretical and experimental investigation of double pass 
finned plate solar air heater was carried out by El-Sebaii et al. 
[34]. The numerical calculations were carried out for Tanta (Lat- 
itude 30° 47’N) prevailing weather conditions. The effect of 
mass flow rates of air on pressure drop, thermal and thermo 
hydraulic efficiencies of double pass finned plate and v cor- 
rugated absorber plate were studied. It was observed that the 
outlet temperature of DPVCPSAH is 2.1-9.7% higher than that 
of DPFIPSAH. DPVCPSAH was 9.3-11.9% efficient than that of 
DPFIPSAH and maximum thermo hydraulic efficiency of DPVCP- 
SAH was 17.4% higher than that of DPFIPSAH and optimum 
values of thermo hydraulic efficiency for both the system are 
obtained at 0.0125 and 0.0225 kg/s. El-Sebaii et al. [35] carried out 
experimental and theoretical investigation of double pass flat and 
v-corrugated plate solar air heater for Tanta (Latitude 30° 47’N) pre- 
vailing weather conditions. The effect of mass flow rate of air on 
pressure drop, thermal and thermo hydraulic efficiencies of dou- 
ble pass flat plate and v corrugated absorber plate was studied. 
Study reveals that the outlet temperature of DPVCPSAH is 5% higher 
than that of DPFPSAH and DPVCPSAH is 11-14% more efficient than 
that of DPFPSAH. The thermo hydraulic efficiency of DPVCPSAH was 
14% higher than that of DPFPSAH and optimum values of thermo 
hydraulic efficiency for both the system are obtained at 0.02 kg/s. 


7. Conclusion 


Based on the review of the literature on double pass solar air 
heater system, it has been found that double pass solar air heater 
were widely investigated both analytically and experimentally. A 
number of studies have been carried out in order to investigate the 
effect of various parameters on the performance of double pass 
solar air heater. The mass flow rate of air and packing material 
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Fig. 10. Double pass channel solar air heater. Three arrangements used for absorber plate (I) cans with staggered arrangement (II) cans with inline arrangement (III) plane 


duct without cans [33]. 
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porosity are considered the important parameters that affect the 
performance of the double pass solar air heater packed bed sys- 
tems. The recycle ratio is found an important parameter that affect 
the performance of double pass solar air heater with external recy- 
cle and it observed that the maximum efficiency is obtained if the 
channel depths and mass flow rate of air is same in both the upper 
and lower duct. A number of analytical and experimental stud- 
ies have been carried out with packed bed, fins integrated double 
pass solar air heater which shows significant increase of the per- 
formance compared to the conventional system. Few studies have 
been reported with corrugated absorber surface. Further, no study 
with artificially roughened double pass solar air heater is reported 
so far in order to compare the heat transfer characteristics and fric- 
tion factor with conventional system. There is tremendous scope for 
future study of double pass solar air heater integrated with rough- 
ness element on both upper and lower side of the absorber surface. 
The information presented here will be beneficial for beginners in 
this area of research. 
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